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ABSTRACT 

We present projected constraints on the cosmic string tension, G/i/c 2 , that could be achieved by fu- 
ture gravitational wave detection experiments and express our results as semi-analytic relations of the 
form Gfj,(fl gw h 2 )/c 2 , to allow for direct computation of the tension constraints for future experiments. 
These results can be applied to new constraints on f2 gw /i 2 as they are imposed. Experiments operating 
in different frequency bands probe different parts of the gravitational wave spectrum of a cosmic string 
network and are sensitive to different uncertainties in the underlying cosmic string model parameters. 
We compute the gravitational wave spectra of cosmic string networks based on the one-scale model, 
covering all the parameter space accessed by each experiment which is strongly dependent on the birth 
scale of loops relative to the horizon, a. The upper limits on the string tension avoid any assumptions 
on the model parameters. We perform this investigation for Pulsar Timing Array experiments of 
different durations as well as ground-based and space-borne interferometric detectors. 
Subject headings: early universe — gravitational waves — instrumentation: miscellaneous — methods: 
numerical 



1. INTRODUCTION 

Cosmic strings are one-dimensional topological defects 
of cosmological size, expected to form either during sym- 
metry breaking phase transitions in the early Universe 
(field theory strings) (jKibblel 1976) or at the end of 
inflation in brane-wq rld scenarios (cosmic superstrings) 
(|Sarangi fc Tvdl200l) . 

The energy scale of cosmic strings is defined by their 
linear energy density or tension, (i, usually referenced 
through the dimensionless quantity G^jl/c 2 , where G is 
Newton's constant and c the speed of light. In the case 
of field theory strings, the tension is directly related to 
the energy scale rj of the phase transition that created 
them, fi ~ t] 2 in natural units, % = c = k = 1. In 
the case of cosmic superstrings, their tension is directly 
related to the fundamental string coupling g s and the 
compactification or warping scales. Therefore, the de- 
tection, or at least the constraining of the cosmic string 
tension, provides a unique laboratory for physics at high 
energies. 

A cosmic string network consists of infinite strings 
which stretch beyond ou r horizon and string loops 
(|Vilenkin k. Shellardlll994[) . After its creation, the net- 
work evolves along with the expansion of the Universe 
and is expected to settle into a scaling regime where all 
its fundamental properties scale along with the horizon 
radius ~ t, which is called the one-scale model. The 
scaling evolution of the network is attained through the 
creation of loops. When two cosmic strings intersect, 
they intercommute with a probability p and exchange 
partners, chopping off loo ps from the network (|Shellardl 
ll987tlJackson et al1l2005h . In the one-scale model, loops 
are born with a size £b equal to a fraction of the parti- 
cle horizon dn at the time of their birth, £b(t) — adn(t), 
where a is a constant. These loops subsequently decay by 
emitting all of their energy in various forms of radiation, 



with the most domin ant of them exp ected to be gravita- 
tional waves (GWs) HVilenkirJ ([1981 : see. lVincent elaTl 
(|1997f ) for an alternative). The GW emission from all the 
loops created by the network creates a stochastic gravi- 
tational wa ve background (SGWB). 

Recently, iSanidas et al.1 (|2012f ) have studied in detail 
the properties of the cosmic string SGWB and set a 
conservative upper limit on the cosmic string tension, 
G/i/c 2 < 5.3 x 1CP 7 for networks with a > 10~ 9 , avoid- 
ing any assumption on the parameters describing the 
GW emission mechanism. Our result is approximately 
one order of magnitude weaker than the e quivalent limit 
found by the iDamour &: Vilenkinl (|2005l ) analytic rela- 
tion Gfj,/c 2 < 1.2 x 1CT 8 (jvan Haasteren et al.l [20TTL 
l2012t ). which is due to the stron g assumptions mad e for 
the latter (see also discussion in lShlaer et al.l (|2012D ). 

As part of the comm on effort of the Europe an Pulsar 
Timing Array (EPTAFI. IPerdman et all (|2010D ) collabo- 
ration to detect GWs with pulsar timing, we present here 
an expansion of our previous work to evaluate the ten- 
sion constraints that will be achieved by future Pulsar 
Timing Array (PTA) experiments of different durations 
and ground-based/space-borne interferometers. This not 
only allows us to predict future constraints, but the semi- 
analytic formulae we deduce can be used when the lim- 
its on f2 gw /i 2 become available. First, we establish the 
parameter space accessible by each of these detectors. 
Then, we discuss how the exclusion curves on the pa- 
rameter space evolve relative to Q gw h 2 . Using this in- 
formation, we will evaluate the tension constraints as a 
function of fl gw h 2 for a range of experiments. We will 
also derive semi-analytic formulae for ease of use in future 
investigations. Additionally, we discuss how the radiative 
efficiency parameter T, affects these constraints. 

1 http://www.epta.eu.org 
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2. CONSTRAINING THE STRING TENSION 

2.1. The cosmic string SGWB 

In the one-scale model, the basic parameters that de- 
termine the GW spectrum are: the cosmic string tension, 
G/i/c 2 ; the birth scale of loops relative to the horizon, 
a; and the intercommutation probability, p. In order to 
compute the SGWB, we also need to model the spec- 
trum of radiation emitted by a cosmic string loop. This 
is quantified by: its spectr al index q, with q = 4/3 for 
cusps and q = 2 for kinks (|Vilenkin fc Shellardlll994T ): a 
cut-off on the number of emission modes, = 1 — > oo, 
which i s used to model the e ffects of gravitational backre- 
action ijCaldwell et al.l ll 9 96) ; and the radiative efficiency 
r « 50, which d escribes the efficiency of the GW emis- 
sion m echanism ()Casper fc Allenlll995l) . In lSanidas et al.l 
( 2012) we have described in detail the computation of the 
cosmic string SGWB which we will only briefly recap it 
here. 

From the energy-momentum tensor for cosmic strings 
wc find the energy lost into loops in order for the network 
to maintain scaling, 



parameter at the present time to, is then given by 
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(1) 

where Poo(i) = A/id^ 2 (t)c 2 is the energy density of in- 
finite strings, A is the number of infinite strings within 
the horizon, (v 2 ) is the mean squared velocity of cosmic 
strings, a(t) the scale factor and V = a 3 (t) is the compu- 
tation volume. We have considered A and (v 2 ) as fixed 
parameters, using their mo st recently published values 
(|Blanco-Pillado et al.ll201lD . quoted in Table [TJ Assum- 
ing that loops have a length i\,{t) = adn(t) when they 
are born, from Eq. ([I]) we can find the loop formation 
rate, 



loop 
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(2) 



where N\ oov is the total number of loops within V{t) 
created since the creation of the network at tf = 
tp\c 4 /(G[i) 2 , with tpi the Planck time. Since the rate of 
decay of length for loops is TGli/c, we are able to com- 
pute the number density n(£, t)d£ of loops with length 
between I and I + di present at time t. 

In our GW emission model, each loop emits into a 
series of harmonics n, with frequency 



_ 2nc 
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and power dE gWi \ oop /dt = P n Gii 2 c with P n = 
Fn~ q / J2m = i m ~ q ■ Then, using Eq. ([3]) and carefully 
redshifting the frequencies to the present time we can 
construct %(/, t) that provides the number density of 
loops, which at the jth mode and at time t, emit GWs 
detected at frequency /. The computation of the en- 
ergy density of GWs per logarithmic frequency interval 

n gw h 2 , where H = lOO/ikms^Mpc- 1 is the Hubble will discuss our choice later. 



(4) 

where p cr i t = 3HyC 2 /8irG is the critical energy density. 

A correction has to be applied in Eq. @ due to the an- 
nihilation of massive particles during the radiation era. 
When the temperature of the Universe, T, drops below a 
specific particle mass threshold, the relevant family be- 
comes non-relativistic, changing the number of relativis- 
tic degrees of freedom, g*. This will result in a faster 
expansion of the Universe every time suc h a mass thresh- 
old is crossed, affecting the value of a(t) (jKolb fc Turnerl 
1990). Although these changes will not significantly af- 
fect the evolution of the cosmic string network, which 
will quickly adapt to the new expansion rate and con- 
verge to its scaling evolution, they will have a significant 
effect on the energy density of the GWs emitted by the 
network. 

These effects can be included with the use of a multi- 
plicative factor for p gw , and therefore f2 gw /i 2 , which has 

the form (g*/gi) 1 ^ 3 , where g^ and g* are the relativistic 
degrees of freedom before and after the transit ion, and 
acts to decrease the amplitude of the SGWB ([Bennett] 
fl986t [Caldwell fc Allerifl99"l IBinetruv et al.ll2012D . This 
correction is applied at time t cor = (327rG j o/3)^ 1 / 2 , 
where p = 7r 2 g»T 4 /30 and its effects are observed at 
a frequency / cor w 2a(t COI )c[f r ad H (i C or)a(io)]~\ which 
depends strongly on a. PTAs, probing nHz GWs, are 
not significantly affected by these corrections, and that 
only in the case where a is large. The se corrections were 
not considered in lSanidas et al.l (|2012t ). but this does not 
affect the robustness of the results presented in the afore- 
mentioned paper. However, when considering GW de- 
tectors operating at the Hz and mHz bands significant 
corrections are expected. We have computed these cor- 
rections for every p article in the S t andar d Model using 
the mass values in iBeringer et all (|2012t l and incorpo- 
rated them to the computed spectra. In Fig. Q] we present 
a typical cosmic string GW spectrum for a = 0.1 with 
the characteristic peak and flat part originating from the 
loop populations in the matter and radiation eras respec- 
tively, along with the corrections to the spectrum due to 
particle annihilation. 

2.2. Exclusion Curves 

The computation of the tension exclusion curves and 
upper bounds for a SGWB limit i s don e using the pro- 
cedure described in lSanidas et al.l (|2012t ). First, we com- 
pute the GW spectra of cosmic string models covering all 
the accessible parameter space. Second, we find which 
configurations are compliant with the SGWB limit, both 
in terms of amplitude and assuming a flat local slopeEL 
Finally, we plot these configurations in the Gfi — a pa- 
rameter space and the maximum tension allowed for a 
given J7 gw /i 2 limit will provide the tension upper limit. 

For a particular experiment, it is essential to determine 
the range of a values which are accessible to it. The 

2 Since we are only considering projected constraints by future 
experiments, an assumption for the spectral slope is necessary. We 
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SGWB of cosmic strings can span a wide range of fre- 
quencies, from below the nHz regime to beyond a THz, 
which depends on the size of the loops created by the 
network. The minimum frequency / m ; n of the SGWB is 
provided by the first emission mode of the loops which 
are born at to, and hence a m ; n = / m in/r<fe(io) /2c. 

In this work, we will not treat the GW detectors as 
broadband, but instead we will designate each one a fidu- 
cial frequency; the one at which maximum sensitivity is 
achieved. For ground-based interferometers (hereafter, 
"Hz-detectors" ) and space-borne interferometers (here- 
after, "mHz-detectors"), we have assumed fiducial fre- 
quencies of / = 100 Hz and / = 0.01 Hz respectively. 
For PTA experiments, the fiducial frequency is the in- 
verse of the elapsed time over which observations were 
taken: / = 6.3 nHz, / = 3.2 nHz and / = 1.6 nHz 
for 5-, 10-, and 20-year PTA experiments respectively. 
They can only detect GW emission for a > a m ; n , where 
log 10 a min = -9.5 for a 5-year PTA, log 10 a min = -9.2 
for a 10-year PTA, log 10 a min = -8.9 for a 20-year PTA, 
log 10 a m in = —15.7 for mHz-detectors and log 10 a m ; n = 
— 19.7 for Hz-detectors. Of course, real detectors are 
broadband and can observe at higher frequencies than 
the one that maximum sensitivity is achieved, but, since 
they are typically only sensitive for a few orders of magni- 
tude in frequency, a m ; n provides a useful rule-of-thumb. 

In Fig. [5] we present the exclusion curves for net- 
works with p = 1 for three projected limits on Q gw h 2 at 
/ = 6.3 nHz and / = 100 Hz. In the case of PTAs (upper 
panel), for fl gw h 2 > 10~ 10 , the upper limit on G/i/c 2 is 
provided by the peak in the mid-a region. The rise in 
the exclusion curve for a ~ a nml is expected, since in 
that region the detector probes the spectrum at frequen- 
cies lower than the peak frequency, with higher values of 
G/i/c 2 required for a fixed Cl gw h 2 . For fi gw /i 2 <C 10~ 10 , 
the characteristic peak disappears and the upper limit 
on the tension is provided by the models with a — a m ; n . 

The Qg W h 2 value at the low frequency cut-off is pro- 
vided by the loops which are born at the present time and 
therefore, it will not be stochastic. In reality, the SGWB 
is expected to start from a frequency / > /mm, where the 
individual loop emission becomes unresolvable, reducing 
the height of the secondary peak and overestimating the 
tension constraints for il gw h 2 < 10~ 10 . However, since in 
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Fig. 2. — Exclusion limits for cosmic string networks with p = 1 
and various values for the upper limit on f2 gw /t 2 at / = 6.3 nHz 
(upper panel) and / = 100 Hz (lower panel). 

that case the spectrum has a positive slope, our assump- 
tion for a flat spectrum leads to an underestimation of 
the tension constraints. These two effects are expected 
to be similar and cancel out, preserving the robustness 
of our results. 

The shape of the exclusion curves for Hz- and mHz- 
detectors exhibit a much less prominent mid-a peak (see, 
for example, Fig. [5J lower panel) . The whole mid-a re- 
gion would be flat if we had not considered the mas- 
sive particle annihilation corrections; a result of the fact 
that such detectors probe only the flat part of the spec- 
trum originating from the radiation era. Moreover, the 
upper limits on the tension are always provided by the 
mid-a peak since the peak created by the networks with 
ol ~ a m j n never gets higher. 

3. PROJECTED CONSTRAINTS 

We computed the spectra for ~ 32000 different pa- 
rameter sets (~ 1.3 x 10 5 C PU hours), using an u pdated 
version of the code used in iSanidas et al.l (|2012D . Most 
changes concern numerical accuracy improvements and 
the use of the recently published values for the infinite 
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TABLE 1 

Parameter values 
used/investigated in this work. 



Quantity 


Value(s) 


A (radiation era) 


45 


A (matter era) 


35 


(v 2 )/c 2 (radiation era) 


0.40 


(v 2 )/c 2 (matter era) 


0.35 


ft 


0.71 


Gn/c 2 


[io- 13 ,io- 4 ] 


a 


[io- 197 ,o.i] 


1 


-4/3,-2 


n* 


1,10 4 


r 


50, [10, 100] 



string parameters from iBlanco-Pillado et al.l ()2011[ ) . A 
list of all these parameter values along with the value 
ranges we covered for the fundamental cosmic string 
model parameters is presented in Table [T] It is only nec- 
essary to use a small number of spectral parameters, since 
the tension constraints for any a in the Gfi/c 2 — a param- 
eter space will always be provided by the exclusion curves 
of either the q = 4/ 3, n* = 1 or the q = 4/3, n* = 10 4 
models as shown in iSanidas et a l. (2012). Additionally, 
we only considered networks with p = 1 since the final 
results can be easily rescaled for networks with p ^ 1. 
For the core results of this paper we used r = 50 
(|Casper fc Ailed Il995h . 

We have investigated SGWB amplitudes in the range 
n gw h 2 £ [10- 12 - 5 , 10~ 7 ], with Alog 10 (fi gw /i 2 ) = O.l. For 
each value for Q gw h 2 , we have computed the exclusion 
curves as a function of a and located the highest pos- 
sible G/i/c 2 value. In Fig. [3] (left panel) we present the 
projected tension constraints as a function of tt gw h 2 . 

For PTAs, the constraints on the tension are almost 
independent of the frequency of maximum sensitivity ex- 
cept in the region 10 -10 < f2 gw /i 2 < 10 -9 . Longer du- 
ration experiments are preferred in that region^, which 
is the expected sensitivit y to the SGWB for LEAP 
(|Kramer fc Stapperi 120101) . The change in the slope of 
the curve at n gw h 2 ~ 10 -10 signifies the point at which 
the tension upper limit is no longer provided by the mid- 
a peak but instead from the a ~ a m i n models. The mHz- 
detectors will provide weaker constraints on the tension 
for all fi gw /i 2 values investigated. The same applies for 
Hz-detector, which will provide weaker constraints than 
both PTAs and mHz-experiments. 

The relatively simple origin of the constraint curves 
makes it possible to describe them with fitting formulae 
Gfi(il gw h 2 ) /c 2 . This will allow any future constraint on 
Q gw h 2 by any experiment operating in the relevant fre- 
quency range to be directly connected to a robust con- 
straint on the string tension. Using least squares fitting 
and setting x = log 10 (f2 gw /i 2 ), we find that: 

3 Longer duration experiments will eventually collect more data 
and set more stringent constraints on f2 gw /i 2 than the shorter du- 
ration ones. However, if two such experiments manage to set the 
same constraint on O gw /i 2 , the longer duration one will provide 
better constraints on the string tension in that region. The bet- 
ter performance of longer duration experiments stems from the 
lower frequency that they probe with maximum sensitivity, plac- 
ing them "deeper" in the high amplitude, matter era peak of the 
cosmic string GW spectrum. 



For 5-year PTA experiments, 

f 0.0959cc 3 + 2.292x 2 + 19.23a; + 50.31 
for - 10 < x < -7 , 



log 10 (G/i/c 2 ) < ^ 



0.0917a; 2 + 2.85a; + 10.5 
Uor - 12.5 < x < -10. 



For 10-year PTA experiments, 

f 0.126a; 3 + 2.972a; 2 + 24.36a; + 63.22 
for -9.8 < x < -7, 



(5) 



log 10 (G M /c 2 ) 



0.0794a; 2 + 2.53a; + 8.52 
Uor - 12.5 < x < -9.8. 



(6) 



For 20-ycar PTA experiments, 

f 0. 1565a; 3 + 3.629a; 2 + 29.05a; + 74.31 
for - 9.5 < x < -7, 



logio(GWc 2 ) = < 



0.074a; 2 + 2.37a; + 7.46 
Uor - 12.5 < x < -9.5. 



For mHz-detectors, 



(7) 



log 10 (G/x/c 2 ) = 1.09a; + 3.16 for - 12.5 < x < -7 . (8) 
For Hz-detectors, 

\og w (Gn/c 2 ) = 1.03a; + 2.83 for - 12.5 < x < -7. (9) 

These formulae can be easily adapted t o describe cos- 
mic s tring networks with p 1. In Sanidas et al. 
(2012) we have pointed out that a non-zero intercom- 
mutation probability only rescales the GW spectrum, 
fl gv/ h 2 oc p~ k , where k is an index describing how the in- 
tercommutation probability affects the population statis- 
tics of infinite strin gs, a nd whose value is discussed in 
i Sakellariadoul (|2005l ) and lAvgoustidis fc Shellardl (|2005l 
2006). Therefore, Eqs. ©, @, ©, © and © can pro- 
vide constraints for any possible combination of p simply 
by multiplying the polynomials with p~ k . 

Additionally, we have investigated the effects of vary- 
ing r on the tension constraints. The value of T depends 
on the trajectory of each loop and is expected to be in the 
range ~ 50 - 100 (lAlTen fc CasperllT99l ICasper fc Ailed 
fl99a iBohg [20111: ICopi fc Vachaspatil 1201 ID . We have 
computed the tension constraints for networks with V € 
[10, 100] and we present some typical results in the right 
panel of Fig. [3] In the case of Hz- and mHz detectors 
the effects of varying T do not significantly change the 
tension constraints. The same holds true for PTAs, ex- 
cept in the region 10 -10 5 < Vt gvJ h 2 < 10~ 9 where low 
values of P are equivalent to decreasing the frequency of 
maximum sensitivity. 

We will conclude this section by discussing the con- 
sequences of our assumption of a locally flat spectrum. 
For Hz- and mHz-detectors, this assumption holds true 
since these detectors probe the flat part of the spectrum 
for most of the G/i — a parameter space. In the case of 
PTAs, the slope will be negative when £l gw h 2 > 10~ 10 
(upper limit provided by the mid-a peak) and positive 
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Fig. 3. — Cosmic string tension upper limits for networks with p - 
panel) and the effects of varying T (right panel). 

when Q gw /i 2 < 10~ 10 (upper limit provided by the peak 
at a ~ ck m j n ). In the first case, the fiat spectrum assump- 
tion will lead to an underestimation of the constraint. 
For example, the E PTA limit on the string tension of 
ISanidas et al.l (|2012f) will provide a hint on the expected 
discrepancy. Reworking that analysis under the assump- 
tion of a flat spectrum we find G/i/c 2 < 5.6 x 10~ 7 
whereas the limit we acquired using the extra slope in- 
formation was G/i/c 2 < 5.3 x 10 -7 , a 5.6% discrepancy. 
We anticipate that the projected constraints in the whole 
range where ft gw h 2 > 10 -10 will be higher by a similar 
percentage, something that does not affect the robustness 
of our approach since we are interested in conservative 
constraints. In the case where, Q gw /i 2 < 10~ 10 the con- 
straints are expected to be overestimated by a similar 
percentage. However, expect that this will be canceled 
by a similar order underestimation of the constraints in 
this region due to the assumption of the stochastic nature 
of the background at / m ; n , as discussed in Sec. 12.21 

4. DISCUSSION 

In this paper, we have presented fitting formulae for 
the upper limits on the cosmic string tension which can 
be used in future GW detection experiments. These con- 
straints are free of any assumption on model parame- 
ters within the one-scale model. These formulae provide 
a conservative alternative to the approach discussed in 
iDamour fc Vilenkml (|2005| ) which is based on cusp emis- 
sion and an approximate loop number density, assump- 
tions which are much stronger. 

The GW spectum of cosmic strings, being extremely 
broadband is potentially accessed by every GW detec- 
tion experiment, with each of them probing bands with 
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logmll^/! 2 

1 as a function of f2 gw /i 2 for various GW detection experiments (left 

different properties. PTAs, probing the high amplitude 
matter era peak of the spectrum for most cosmic string 
model parameter combinations and only partially suffer- 
ing from the massive particle annihilation corrections in 
the case of large loops, provide an excellent tool to detect 
strings, or at least, set the most stringent constraints on 
their tension. On the other side, ground-based interfer- 
ometers probe a much wider area of the parameter space 
(own ~ 10 -20 ) making them the only tool able to detect 
emission from very small loops, and they are not sensi- 
tive to the, unknown, details of the radiation spectrum 
(q and n»). Space-borne interferometers, operating in 
a frequency range between PTAs and ground-based in- 
terferometers, share advantages and disadvantages from 
both. 

PTAs already set the most stringent constraints on 
the amplitude of the SGWB when compared to LIGO, 
and Advanced LIGO which is expected to become op- 
erational in 2015 will achieve a similar perfor mance 
(|Harrv fc the LIGO Scientific Collaborationll20iol) . The 
further increase in the sensitivity of PTAs expected in 
the near future (LEAP is expected to set upper lim- 
its on fl gw h 2 < 10~ 10 ), will lead to constraints on 
G/^/c 2 < 2 x 1CT 9 , with a SKA0-scale PTA experiment 
improving this by two orders of magnitude. 
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